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Abstract 

By using the effective potential approximation and taking into ac- 
count the dominant top quark and scalar top quark loops, radiative 
corrections to MSSM Higgs potential are computed in the presence of 
the supersymmetric CP-violating phases. It is found that, the lightest 
Higgs scalar remains essentially CP-even as in the CP-invariant the- 
ory whereas the other two scalars are heavy and do not have definite 
CP properties. The supersymmetric CP-violating phases are shown 
to modify significantly the decay rates of the scalars to fermion pairs. 
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1 Introduction 



The Minimal Supersymmetric Standard Model (MSSM) consists of various 
soft supersymmetry breaking parameters as well as the Higgsino mass pa- 
rameter, fi, coming from the superpotential. In general, there is no a priori 
reason for taking all these parameters real, and thus, Yukawa couplings, gaug- 
ino masses, trilinear Higgs-sfermion couplings, Af, Higgs bilinear coupling, 
m\ = fiB, and fi- parameter itself can all be complex. On the other hand, the 
MSSM Lagrangian has two global symmetries U(1)pq (Peccei-Quinn symme- 
try) and U(1)r-pq (an R symmetry) under which all fields and parameters 
are charged. The selection rules for these symmetries limit the combinations 
of dimensionful parameters that can appear in a physical quantity so that 
one has, in fact, only three of these phases physical |0, Without loss of 
generality, these three physical phases can be identified with (1) the phase in 
the Cabibbo-Kobayashi-Maskawa matrix, 5ckm, (2) = Arg(fi), and (3) 
ip A = Arg(Af). Thus any physical quantity T has an explicit dependence 
on these phases: T = F{&ckMi W V ? a / )- 

Despite their presence in the Lagrangian, the phenomenological relevance 
of these CP-violating phases has often been questioned due to the smallness 
of the neutron and electron electric dipole moments |3| which require 
them to be at most O(10~ 3 ). However, recent studies have shown that it is 
possible to suppress neutron and electron dipole moments without requiring 
these CP-violating phases to be small by allowing the existence of either non- 
universal soft breaking parameters at the unification scale |4j or some kind 
of cancellation among various supersymmetric contributions ||, or heavy 
enough sfermions for the first two generations ||. In fact, following the 
last scenario, it was recently shown that the CP-violation in B- and K- 
systems can be saturated with 92 M and ipA f only [[/J. Apart from electric 
dipole moments and weak decays, these phases play a crucial role in the 
creation of the baryon asymmetry of the universe at the electroweak phase 
transition [Q. 

In this work, assuming that the electric dipole moments are suppressed 
by one of the methods metioned above, we take supersymmetric phases un- 
constrained, and investigate their effects on the Higgs sector of the MSSM. 
At the tree-level the Higgs sector of the MSSM conserves CP due to the fact 
that the superpotential is holomorphic in superfields entailing the absence 
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of flavour changing neutral currents and scalar-pseudoscalar mixings. When 
the supersymmetric phases ip^ and ip^. vanish Higgs sector conserves CP at 
any loop order. In fact, the CP-conserving Higgs sector has been analyzed 
by several authors with the main purpose of evaluating the mass of the light- 
est Higgs boson which has the tree- level upper bound of Mz- It has been 
found that radiative corrections, dominated by top and stop loops, elevate 
the tree-level bound significantly ||. These one-loop results @ have been 
improved by utilizing complete one-loop on-shell renormalization [|H| , renor- 
malization group methods JTT1, diagrammatic methods with leading order 



QCD corrections [12|, and two-loop on-shell renormalization [13]. However, 



when the supersymmetric phases are non-vanishing, as the recent studies 
have shown [[TJj], the Higgs sector become CP- violating through the radia- 
tive corrections. As in the CP-conserving case, the radiative corrections will 
be dominated by the top and stop loops. Below we investigate effects of 
the supersymmetric phases on the Higgs masses, scalar-pseudoscalar mix- 
ings, and decay properties of scalars to fermion pairs. In doing this, we shall 
calculate one-loop radiative corrections coming from top and stop loops in 
the effective potential approximation. 

This work is organized as follows. In Sec. 2 we compute the one-loop 
effective potential using top quark and stop contributions together with the 
specification of the particle spectrum and mixings. In Sec. 3 we discuss, as 
an example, the decay properties of the Higgs scalars to fermion pairs. In 
Sec. 4 we conclude the work. 



2 Effective potential 

The Higgs sector of the MSSM consists of two SU(2) doublets H 1 , H 2 with 
opposite hypercharges Y\ = —1, Y 2 = +1, and non- vanishing vacuum expec- 
tation values Vi, v 2 - Allowing a finite alignment, 6, between the two Higgs 
doublets, we adopt the following decomposition: 

Hi = 
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At the tree level the Higgs sector is described by the scalar potential 



V (H 1} H 2 ) 
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with the parameters 



m i = m \ + > m 2 = + l/^l 2 > m s = l/^l > A i = A 2 = (#2 + #i)/ 4 

,2 „2\ /j \ _ „2 



A12 = (gi-gi)/*, Xn = -9 2 2 /2 (3) 

where m 2 ^ and i? are the soft supersymmetry breaking parameters. As 
is seen from (3), the tree level potential is described by real parameters; 
thus, the alignment between the two doublets can, in fact, be rotated away. 
Since, in the minimum, the potential is to have vanishing gradients in all 
directions, in particular, dVo/d(pi )2 = mgsin^ = 0, one automatically gets 
9 = 0. Those terms of the tree level potential (2) quadratic in the components 
of the Higgs doublets (1) give the mass-squared matrix of neutral scalars the 
diagonalization of which yields the CP=-1 boson A = cosj3ifi — sin/?<^ 2 
with mass M\ Q = — Tn 2 ,/ sin/? cos/?, and two CP=+1 bosons which are linear 
combinations of <p\ and </> 2 with a mixing angle a. The mixing angle a and 
the masses of the CP even scalars h and H are given by |L5| 



M% +Ml 

tan 2a = — f -f tan 20 4 

M y - Mf 

1 



M 2 h{H) = 2^ + M z -{+)\j{M 2 AQ + Miy-AMl Ml cos 2 2(3) (5) 

where tan/? = v 2 /vi, and Mf = {g\ + g\){y\ + f|)/4 in our convention. It is 
readily seen that for tan/? > 2 one has /? ~ tt/2, M h ~ M z and M# ~ M^. 
However, it is known that radiative corrections elevate (bounded by Mz 
at the tree level) significantly || without modifying the mass degeneracy 
between H and A when the theory conserves CP. When, however, the CP- 
violating MSSM phases are switched on, the degeneracy between H and A 



can be lifted considerably as discussed in [£L4[ . 

We now start computing the radiative corrections to the tree potential 
(2) in the presence of the CP-violating MSSM phases f^ t A r Our main 
concern will be the investigation of the masses and mixings of the scalars as 
a function of the CP-violating angles. These mixings could be CP-conserving 
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(like h — H mixing in the CP-respecting limit) as well as CP-violating as 
we will discuss below. To evaluate the radiative corrections we follow the 
effective potential approximation where the tree-level potential (2) is added 
to the one-loop contributions having the famous Coleman- Weinberg |9|, |I~6 
form 

ZV = ^S t rMHHuH 2W o g ^^l>-l) (6) 

where Str = Y,j(~ 1) 2J+1 (2J + 1) is the usual supertrace, and JM.(H\, H 2 ) 
is the Higgs field dependent mass matrix of particles. AV depends on the 
renormalization scale Q which is presumably around the weak scale. In 
evaluating AV one includes the contributions of vector bosons, Higg bosons 
and fermions as well as their supersymmetric partners gauginos, Higgsinos 
and sfermions. Among all these particles top quarks and scalar top quarks 
give the dominant contributions |§. However, for very large tan/3 values 
bottom-sbottom and r-stau systems can become important. Besides these, 
since the dependence of on the CP-violating phases f^Af originates 
from only the Higgsino (through fi dependence) and sfermion (through /i and 
Af dependence) mass matrices, these two particle species attain a seperate 
importance. However, if one wishes to include the Higgsino contribution, all 
particle species must be included since then precision of the computation rises 
to the level of gauge couplings. In the following we neglect the contributions 
of gauge couplings and restrict ourselves to moderate values of tan/3 so that, 
to a good approximation, the dominant terms in AV^ are given by top- stop 
system. This approximation is convenient in that it picks up the phase- 
sensitive dominant contributions to AV^. 

In (tL,tu) basis stop mass-squared matrix, neglecting the D-term contri- 
butions, takes the form 

, M| + h*\H$\ 2 h t {AtH%-ii*Hr, i m 
' 1 h t (AtH? - /ifl?) MI + /M 2 ; [<) 

where /i and A t are complex, M| ^ are the soft mass-squareds of left- and 
right-handed stops, and ht is the top Yukawa coupling. Denoting the eigen- 
values of Mi by m| and using m\ = h 2 t \H$\ 2 the one-loop effective potential 
takes the form 

V = V, + JL, ( £ ro >g f 2 - \) - 2 m? (.og f - §)) . (8) 

\a=ti,t2 I 
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We require this effective potential to be minimized at (vi,v 2 ,6) at which it 
has to have vanishing gradients in all directions and the masses of the Higgs 
scalars must be real positive. Gradients of the potential with respect to the 
charged components of the Higgs doublets automatically vanish as there is 
no charge breaking effects in the vacuum. On the other hand, extremization 
of V with respect to neutral components of the Higgs doublets yield 

ui {2m\ + Ai^i + (Ai 2 + Ai 2 )u|) + 2m\v 2 cos 6 + 2 = (9) 

v 2 (2ml + X 2 vl + (X 12 + X 12 )vl) + 2mlv lC os6 + 2\^-] = (10) 

\ 0(f) 2 J 

mj^i sin - ^^-^ =0 (12) 

where the subscript "0" implies the substitution (pi = <p 2 = ipi = ip 2 = in 
the corresponding quantity. Equations (9) and (10) come by no surprise as 
they are the counterparts of the ones occurring in the CP-conserving case. 
However, with (dAV/dipi t2 ) Q ^ 0, equations (11) and (12) now imply a non- 
trivial solution for sin 6 unlike the CP-conserving case where these gradients 
indentically vanish and one automatically obtains a vanishing 8. After some 
algebra one can show that 

fdAV\ n fdAV\ 

Therefore, equations (11) and (12) imply one and the same solution for 9 

m 2 3 sin6 = ^f3 ht \fi\\A t \ sin 7/(7^, m|) (14) 
where (3 ht = (3h 2 t ) / 'lQn 2 , 7 = <p„ + (p At , and 

f( \ o , 1 xy , y + x * y f,*\ 

f(x, y) = -2 + log — + log- 15 

Q A y - x x 

is a scale-dependent one-loop function. Actually, if one uses the decompo- 
sition of the Higgs doublets in (1), in all one-loop formulae 7 gets replaced 
by 7 + 9. However, since 6 itself is a loop-induced quantity, its appearence 
together with 7 is a two- and higher- loop effect which we neglect. Thus, 
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when computing Ay we drop 9 from (1) knowing that it is induced through 
(14). It is readily seen from (14) that unless 7 vanishes 9 remains finite. 
Furthermore, due to the form of the stop mass-squared matrix, all one-loop 
quantities turn out to depend on the combination 7 = <y? M + (pA t - Of course, 
had we included the Higgsino contributions there would be terms that de- 
pend solely on (p^ destructing this kind of relation. However, they would be 
subleading compared to top and stop contributions discussed here. 

In (14) and all formulae below m| 2 denote stop mass-squared eigenvalues 
evaluated at the minimum of the potential: 

< 2 = ^( M ! + M ! + 2 ™t 2 TA?) (16) 

where 

A- = v /(M|-M|) 2 + 4mf(|A| 2 + |/i| 2 cot 2 /^-2|/i||A|cot/5cos7) . (17) 

As usual, construction of the mass-squared matrix of the Higgs scalars 
proceeds through the evaluation of 

/ d 2 V \ 

M2= \dx~dx~) ' whereXi E B = ifa'fc'Vi'V*} ■ ( 18 ) 

Solving m\ 2 from the stationarity conditions (9) and (10), and replacing 
them in M 2 one observes that, in the basis £>, the vector {0, 0, — cos (3, sin (3} 
corresponds to the Goldstone mode G° eaten by Z boson to acquire its mass. 
Then in the reduced basis £>' = {0i, 2 , sin (3(fi + cos /S^} the mass-squared 
matrix of the Higgs scalars becomes 

/ M 2 z c 2 +M 2 A s} + A n -(M 2 +M 2 A ) S/3 cp + A 12 rA \ 
M 2 =\ -(Mf + M 2 A )s p cp + A 12 M 2 s 2 + M\c 2 + A 22 sA (19) 
V rA sA M\ + A J 

where cp = cos/3 and sp = sin/3. The scalar mass-squared matrix involves 
various parameters whose explicit expressions we list below. The adimen- 
sional parameters r and s are given by 



sin f3\A t \ 


cos 7 — 


cot (3 


sin 7 


M 





Sill" 

sin (3 



sin7 \/i\\A t \g(ml,m 2 i2 ) 



(\A t \(\A t \ - \fi\ cot (3 cos7)g(m h ,m t2 , 
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(21) 
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While r originates mainly from the stop left-right mixings, s has an additional 
term depending on the the stop mass splitting log — 

The parameters of mass dimension can be expressed as follows 



sin p [m~ — m z 
A O o \f^\ 2 m 2 t {\A t \ cos-f - \fi\ cot (3f , ., 



m~ — m~ 

12 tl 



A no ii 2 rl^l cos 7- cot/3 m| 
A12 = -2/3 ht |//|m t { 7—= log - - 



t 2 t\> tl 



\A t \[(\A t \ cos7 - |/x| cot /?) 2 + IAKIAI - Mcot/3) sin 2 7] 

(m| -m|) 2 

x ^K,m|,)} (24) 



A 2.-? m 2 rinr m ^ I 2 l^|(|A|-|/i|cot/3cQS7) 
A22 - ^Ph t m t {log t 1 -, — 5 5\" lo & 



n? 



*2 



4 1 I 2 2 \ ivy fc> 2 

|AP(|^|-Mco^oos 7 )^ K ; m|)} ' ' (25) 



v *2 tl' 



where the function gr(m| , m| ) in these expressions reads 

tffoy) = /v>,y) -loggj- (26) 

Therefore, unlike /(m^m- ), g{m~ ,m~) does not have an explicit depen- 
dence on the renormalization scale Q- In fact, the adimensional parameters 
r and s as well as the mass parameters (21)-(24) have no explicit dependence 
on Q. On the other hand, the remaining mass parameter, M\, in the scalar 
mass-squared matrix (19) is an explicit function of Q: 

M\ = . - 3 - — 1 . 27 

sin p cos p sin 7 

It is the 9 (13) dependence of M\ that makes it Q-dependent. However, the 
explicit Q dependence of 6 should cancel with the implicit Q dependence of 
tan p and m 2 to make M\ scale-independent || . 

According to the decomposition of Higgs doublets in (1), (f>i and 02 are 
of CP=+1 whereas sin ftipi + cos /3y2 2 is of CP=-1. As suggested by the form 
of the scalar mass-squared matrix (19) there are mainly two kinds of mix- 
ings: (1) mixing of the scalars with different CP induced by Mf 3 = rA 
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and M| 3 = sA, and (2) mixing of the CP=+1 scalars through M\ 2 = 
— (Mf + M'fyspCp + A 12 . While the former are induced purely by the non- 
vanishing supersymmetric phases the latter exists in the CP-respecting limit 
too. In the CP-conserving limit, that is, sin 7 — > 0, one obtains rA — > 0, 
sA — > and A — > 0. In this case CP=+1 and CP=-1 sectors in (19) decouple 
and reproduce the particle spectrum of the CP- conserving limit in which 
Ma becomes the radiatively corrected pseudoscalar mass, and A 11j12j2 2 be- 
come the usual one-loop contributions || to the CP=+1 scalar mass-squared 
matrix. For a proper interpretation of results of the numerical analysis be- 
low it is convenient to know the relative strengths of the CP-violating and 
conserving mixings. For large tan/? equations (20) and (21) go over to 

1 \A t \ 4K + M2) 

r~— cot 7, s~ (-p-p H 1 I, . 1 ) (28) 

SU17 \n\ \^\\A t \ 

where we assumed M^ ~ M^ = Mq » m t \A t \ and \A t \ » |/x| cot /? in 
derivation. Equation (28), thus, implies that |r|/|s| << 1. This follows 
mainly from the dependence of the stop mass-matrix (7) on which causes 
not only |sA| but also IA12I and IA22I to be larger than |An| and |rA| 
through stop and stop-top splittings. An immediate consequence of (28) is 
that one eigenstate of the scalar mass-squared matrix (19) will be of mainly 
CP=+1. Thus one expects that among the three mass eigenstate scalars 
one will continue to have CP=+1 with a small CP=-1 component while the 
other two can mix significantly depending on the relative stregths of the 
other one-loop corrections. This observation can be justified numerically by 
analyzing the relative strengths of CP-violating and CP-conserving mixings 
as classified above. Fig.l shows the variation of Mf 3 /\Mf 2 \ (solid curve) 
and Mf 3 /|M? 2 | (dashed curve) with tan/3 for M z = M^ = \A t \ = 10 • M z , 
\fi\ = 2.5 • Mz and Ma = 2 • Mz with 7 = 7r/4. As the figure suggests larger 
the tan (3 bigger the CP-violating mixings compared to the mixing between 
the CP=+1 components. The increase of these ratios with tan/3 can be 
understood as follows: (1) (3 — > n/2 as tan/5 grows to higher values, and 
thus, I — (Mf + M\)s/3Cp\ decreases gradually, (2) A12 is positive for these 
parameters (and even for higher values of due to cot/3 suppression) and 
grows with tan/3 due to log m~ /m~ so that Mf 2 decreases with increasing 
tan/?, and the ratios increase gradually since Mf 3 and M| 3 decrease with 
tan/3 more slowly. 
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Figure 1: Variation of Mf 3 /|Mf 2 | (solid curve) and M| 3 /|Mf 2 | (dashed curve) with 
tan (3 for M~ L = M k = \A t \ = 10- Mz, |m| = 2.5-M z and M A = 2-M z with 7 = vr/4. 
The CP- violating mixings become important for large tan (3. 

In Fig. 1 we plot ratios of the CP-violating mixings to CP-conserving 
ones; however, if one plots CP-violating mixings directly, for example in 
units of (10 • M z ) 2 for = 10 • M z , 7 = tt/2 and M A = 5 ■ M z , in absolute 
magnitude, M 2 3 (M| 3 ) starts with ~ 4 ■ 10 -2 (~ 10~ 2 ) at tan/3 = 2 and falls 
down 3 • 10~ 5 (2 • 10 -4 ) at tan/3 ~ 30. These results generally agree with 



those of [0 though the computational schemes are different. 

Since the parameter space is too wide to cover fully, in the following we 
restrict ourselves to the following set 

M L = M k = 500 GeV , | A t \ = 1 TeV , = 250 GeV , M A = 200 GeV (29) 

and vary 7 over its full range. Each time we consider low and high tan/3 
regimes seperately by taking tan/3 = 4 and 30. As tan/3 increases cot/3 
decreases, and this enhances the contribution of the radiative corrections. 
However, variation with tan /3 is not the whole story because even for cot (3 ~> 
0, A, An, A12 are proportional to |/x||A t | so that choice for the latter affects 
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the strength of the radiative corrections. Especially for large tan (3, stop 
masses weakly depend on therefore, these elements of the mass-squared 
matrix become more sensitive to the choice for The parameter set (29) is 
a moderate choice in that it enhances the radiative corrections through large 
\A t \ term without causing too big splittings among the the A coefficients 
in (22)-(25) thanks to the relatively small \fj,\ term. Dependence on the 
parameter M4 is as in the CP-invariant theory, namely, heavy scalars have 
masses around M A . 




Figure 2: Percentage composition of Hi as a function of 7 for tan (5 = 4 (left panel) 
and tan/3 = 30 (right panel). Here (f)±, <p2 and sin/3</?i + cos f3(f2 contributions are 
|7£n| 2 (solid curve), |7£i2| 2 (dashed curve) and |7£i3| 2 (short-dashed curve), in 
percents. Values of the parameters are given in (29). 

In principle one can diagonalize analytically the scalar mass-squared ma- 
trix (19); however, the results will be too complicated to be suggestive. In- 
stead of using such oblique expressions we shall fix the notation for diagonal- 
ization and numerically analyze the results. The scalar mass-squared matrix 
(19) can be diagonalized by a similarity transformation 

TL ■ M 2 ■ 7Z T = diag. (M 2 Hi , M% , M| ) , where TZ ■ 7Z T = 1 (30) 
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Figure 3: Percentage composition of Hi as a function of 7 for tan (3 = 4 (left panel) 
and tan/3 = 30 (right panel). Here (pi, (pi and sin/?y?i + cos f3(pi contributions are 
1 7?-2i 1 2 (solid curve), \lZn\ 2 (dashed curve) and |7?.23| 2 (short-dashed curve), in 
percents. 

where the mass- eigenstate scalar fields are defined by 

V H 3 J \ sin /3</?i + cos j3(p 2 

7- dependence of the elements of TZ is crucial for determining the CP- impu- 
rity of the mass eigenstate scalars Hi. In analyzing TZ we adopt a convention 
such that in the limit of vanishing sin 7 we let H\ — > h, Hi — > H and H3 — > A, 
that is, i/i is the lightest Higgs. We expect results of the CP-invariant the- 
ory be recovered at the CP-conserving points 7 = 0, ir, 2n except for 7- 
dependence of various parameters. 

Depicted in Fig. 2 is the 7- dependence of Hi composition in percents 
for the parameter set in eq. (29). From the left panel we observe that, on the 
average, Hi has ~ 30% 0i and ~ 70% 2 composition. As is immediate from 
the figure, CP=-1 component of Hi is small, in fact, it never exceeds 0.02% 
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in the entire range of 7. From the right panel, however, we observe that 02 
contribution rises near to 100% line, and correspondingly, <pi contribution 
remains below 2.5%. This result is a consequence of /3 — > it/2 limit remines- 
cent from the CP-invariant theory. In this large tan /3 limit, sin (3<pi + cos f3(p 2 
composition of Hi reaches at most to 0.2% over the entire range of 7. It is 
clear that both windows of the figure suggest that the lightest Higgs remains 
essentially a CP-even Higgs scalar for the parameter space in (29). 




Figure 4: Percentage composition of #3 as a function of 7 for tan (3 = 4 (left panel) 
and tan/3 = 30 (right panel). Here <f>i, (p 2 and sin/3<^i + cos (3(p 2 contributions are 
1^-31 1 2 (solid curve), 1 72-32 1 2 (dashed curve) and 1 72.33 1 2 (short-dashed curve), in 
percents. 

In Fig. 3 we show the percentage composition of H 2 as a function of 7 
for tan/3 = 4 (left panel) and tan/3 = 30 (right panel). In agreement with 
the left panel of Fig. 2, for tan/3 = 4 (left panel) H 2 has ~ 70% 0i and 
~ 30% 02 composition. Unlike Hi, however, sin/3</?i + cos/3y?2 composition 
of H 2 becomes as large as 1.3%. As expected, this increase in the CP=-1 
component is compansated by if 3. More spectacular side of Fig. 3 arises for 
large tan (3 (right panel) in which H 2 is seen to gain non-negligible CP-odd 
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Figure 5: Masses of the scalars Hi as a function of 7 for tan/3 = 4 (left panel) and 
tan/3 = 30 (right panel). Here Mjix, Mh 2 and Mh z are shown by solid, dashed, 
and short-dashed curves, respectively. In both panels H\ is the lightest scalar 
whose composition is shown in Fig. 2. 

composition. In accordance with large tan (3 limit described by eq. (28) there 
is a strong dependence 7. Sum of <pi and 2 compositions of H 2 starts from 
7 = at 100% line, and the former is diminished rather fast until 7 = n/2. 
Beyond this point it rises rapidly to 98% line at 7 = ir where its CP=+1 
component completes to 100% in accordance with the CP conservation. One 
notes the complementary behaviour of sin (3ipi + cos (3<p 2 composition which, 
in particular, implies that H 2 is a pure pseudoscalar around 7 = ir/2. 

Depicted in Fig. 4 is the percentage composition of H 3 as a function 
of 7 for tan/3 = 4 (left panel) and tan/3 = 30 (right panel). In agreement 
with the left panels of Figs. 2 and 3, H 3 is almost a pure pseudoscalar for 
tan (3 = 4. On the other hand, for tan (3 = 30 (right panel) H 3 is seen to loose 
its CP-purity in accordance with the right panel of Fig. 3. Thus if 3, except 
for the points discussed above, does not have definite CP characteristics. 
From Fig. 3 and Fig. 4 one concludes that heavy scalars have non-negligible 
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CP-impurity in agreemnet with the results of [Q. 

4>i (solid curve), 4> 2 (dashed curve) and sin/3<£>i + cos/3y?2 (short-dashed 
curve) compositions of the Higgs fields Hi shown in Figs. 2-4 need further 
elaboration. In accordance with Fig. 1 CP-violating mixings are small for 
small tan j3 and one necessarily recovers the results of the CP-invariant the- 
ory where <pi and 4> 2 mix with each other to form h° and H°, and sin/3y?i + 
cos /3if2 is nothing but the psedoscalar A . Thus, in this limit mixings are as 
in the CP-invariant theory as is evident from the left panels of each figure. 
For large tan/3, however, each Higgs field undergoes certain variations in its 
compositions. First of all, dominant 02 composition of the light Higgs is 
easily understandable since the analog of the tree-level Higgs mixing angle 



(4) approaches (3 + n/2 « n for large tan (3 [fT?]- The remaining component 
of Hi comes mainly form (pi since CP-breaking contributions are small for 
this eigenvalue (See (28)). On the other hand, compositions of H 2 and H 3 
are determined from large tan /3 (dominant <pi contribution to their CP-even 
parts) as well as the radiative corrections. As is evident from (28), in large 
tan (3 regime M 13 ~ sin 7 cos 7 and M23 ~ sin 7 so that (for example) the ze- 
roes of the components of H 2 follow certain combinations of these functional 
behaviours. The sharp changes in the 02 and sin /3(pi + cos f3tp 2 compositions 
of H 2 and H 3 at 7 = (for large tan/3) follows from the 7 dependencies 
of M 13 and M 23 (See eq. (28)). Indeed, variation of the compositions near 
7 = behaves roughly as sin 27 ± sin 7 which varies quite fast near the origin. 

Until now in Figs. 2-4 we have discussed the CP properties of the eigen- 
states of the scalar mass-squared matrix (19). Now we analyze the 7 depen- 
dence of the masses of these scalars for identifying their hierarchy. Fig. 5 
shows the 7-dependence of the scalar masses for tan (3 — 4 (left panel) and 
tan/3 = 30 (right panel). It is clear that Hi is the lightest scalar in both 
cases. Moreover, the small gap (at most ~ 20 GeV) between M Ha and M H2 
in the left panel is closed in the right panel where H 2 and H 3 are degenerate 



in mass. This behaviour occurs also in the CP-conserving case [15, pfl due to 



the large value of tan /3. One more thing about Fig. 5 is that the masses of 
if 2 and H 3 are almost completely determined by M4 = 2 • Mz- For instance, 
for M4 = 10 ■ Mz, H 2 and i^3 weigh approximately 10 • Mz- Unlike this 
strong Ma- dependence of Mh s and Mh 2 , mixing among the Higgs scalars 
and lightest Higgs mass depends mainly on tan (3. The 7-dependence of the 
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masses around 7 = 7r (especially in the case of tan /3 — 4) differ from those at 
other CP-conserving points due mainly to the minimization of the light stop 
mass (See eqs. (16) and (17)). Finally, one observes that the light Higgs mass 



is higher than the usual constrained MSSM bounds |18| due to the large value 
of A t and relatively small soft stop masses which is not possible to produce 
through RGE's (except, possibly, with non-universal initial conditions). 

From Figs. 2-5 one concludes that the lightest Higgs scalar remains es- 
sentially a CP=+1 scalar, and the remaining heavy scalars do not possess 
definite CP characteristics. In this sense, MSSM is seen to accomodate a light 
CP=+1 Higgs boson as in the CP-invariant case. In addition, the heavy in- 
definite CP Higgs scalars of the MSSM not only have no correspondent in 
the SM but also are distinguishable from the lightest Higgs due to both their 
masses and CP properties. In the next section we shall discuss decay prop- 
erties of these Higgs scalars to fermions and investigate the deviations from 
the CP-conserving limit. 



3 An example: Decays of the Higgs scalars 
to fermion pairs 

The form of the scalar mass-squared matrix (19) as well as the mixing matrix 
1Z shows clearly the CP-violation in the Higgs sector. In general, all parame- 
ters of the Higgs sector turn out to depend on these CP-violating angles, for 
instance, tan /3, masses of scalars, stop masses are explicit functions of 7. In 
this sense couplings of the Higgs scalars to the MSSM particle spectrum are 
necessarily modified. Therefore, the CP-violating angle 7 can show up, for 
example, in the rates for various collision processes testable at future collid- 
ers. As an example, one can consider an e + e~ collider with sufficiently large 
center-of-mass energy. In such a collider one of the main Higgs production 
mechanisms is the Bjorken process Z* — > ZHi. When the center-of-mass en- 
ergy is varied over a range of values including the masses of the scalars it is 
expected that the cross section, as a function of the invariant mass flow into 
Hi branch, should show three distinct peaks situated at the scalar masses 
Mff r Needless to say, in the CP- conserving case there would be two peaks 
instead of three. 

We now discuss couplings of the scalars to fermion pairs in detail. That 
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the scalars Hi are devoid of definite CP properties influence their couplings 
to fermions significantly. Rephasing the fermion fields appropriately one can 
always make fermion masses real after which coupling of the scalar Hi to 
■u-type quarks, d-type quarks and charged leptons take the form 

H t uu : (V2G F )^ 2 ^-(n t2 + icospn z3l5 ) (32) 
sin id 

Hdd : (v^GV) 1/2 ^(fta+*sin/37^375) (33) 

Hll : (v^GV) 1 / 2 ^-^! +i sin pHa lb ) (34) 
cos p 

where one observes that each coupling picks up a 75 piece showing its CP=-1 
content. Moreover, as the Feynman rules above dictate 75 piece is enhanced 
for large tan (3 (cot f3) for -u-type quarks (d-type quarks and charged leptons). 
To have a quantitative understanding of the effects of 7 on the Higgs decays 
to fermion pairs it is convenient to compute the ratio 

p mi* - ffn ± 0) ( . 

~ fS^R) (35) 

(K iq ) 2 (l - Am}/M" H f/ 2 + a 2 (7^ 3 ) 2 (l - 4m 2 /M 2 J 1 / 2 
(n iq (0)f(l - 4m 2 /M^(0))^ + a 2 (^ 3 (0)) 2 (l - 4m 2 /M^(0))V2 

where q = 2(1) and a/ = cos /3(sin j3) for -u-type quarks (d-type quarks and 
charged leptons). The argument "0" of the quantities in the denominator 
implies the replacement 7 = 0. It is clear that T(Hi — > //) consists of 
phase space factors pertinent to both CP=+1 and CP=-1 cases seperately. 
According to the conventions we apply, for % — 1, 2 {% — 3) only the first 
(second) term survives in the denominator. The CP-violating MSSM phase, 
7 not only functions in creating the additional terms in the Feynman rules 
(31-33) but also affects the couplings and masses themselves. 

We now perform a numerical computation of Rif for the parameter space 
used up to now. In particular, we concentrate on two cases f — b and f — c, 
that is, we consider b- and c- quarks in the computation. Such an analysis will 
be exhaustive as it covers the cases listed (32)- (34). We start by listing the 
quantities (^(O)) 2 necessary for computing Rif in Tables 1 and 2. Below, 
when speaking about the elements of the mixing matrices for 7 = we shall 
always refer to these tables. 

Fig. 6 shows the 7-dependencies of R±b and R\ c for tan/3 = 4 (solid 
curve) and tan/3 = 30 (dashed curve). Let us first discuss Ru (left panel). 
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Figure 6: Ru (left panel) and R\ c (right panel) denned in (34) as a function of 7 
for tan (3 = 4 (solid curve) and tan (3 = 30 (dashed curve) . 

As Fig. 2 shows 0i composition of iJi starts with 0.31 at 7 = (See Table 
1) and falls to 0.24 at 7 = tt. Therefore, R lb falls to 0.24/0.31 = 0.77 around 
7 = 7r. On the other hand, for tan/5 = 30 .Ru, obtains a relatively fast 
variation of with 7. This behaviour of follows from its 0i component 
in the right panel of Fig. 2, which takes the value of 1.1% around 7 = n. 
Therefore, Ru, rises to 0.01/0.02 = 0.5 at 7 = it as follows from Table 2. 
On the other hand, from the right panel of the figure one observes that, R\ c 
remains around its conterpart in the CP-conserving limit. Using Tables 1-2, 
and 02 compositions in Fig. 2 one can easily infer the behaviour of R\ c . In 
particular, constancy of tan/5 = 30 curve follows from the constancy of the 
02 composition in Fig. 2 (right panel). 

Depicted in Fig. 7 is the dependence of R 2 b (left panel) and R 2c (right 
panel) on 7 for tan/3 = 4 (solid curve) and tan/3 = 30 (dashed curve). First 
concentrating on R 2 b, one observes a slow variation with 7 compared to the 
lightest Higgs (Fig. 6 (left panel)). The behaviour of tan/3 = 4 curve simply 
follows from the 0i composition of H 2 in Fig. 3 (left panel). However, the flat 
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(^•(o)) 2 


i=l 


i=2 


i=3 


(^i(O)) 2 


0.31 


0.69 


0.0 


(^(o)) 2 


0.69 


0.31 


0.0 


(^a(O)) 2 


0.0 


0.0 


1.0 



Table 1: Elements of the scalar mixing matrix entering the Feynman rules (32)-(34) 
for tan (3 = 4 and 7=0. Here i (j) runs over Hi (0i, 02, sin/3^i + cos (3ip 2 ). 



(^•(o)) 2 


i=l 


i=2 


i=3 


(^i(O)) 2 


0.02 


0.98 


0.0 


(^ 2 (0)) 2 


0.98 


0.02 


0.0 


(^3(0)) 2 


0.0 


0.0 


1.0 



Table 2: Elements of the scalar mixing matrix entering the Feynman rules (32)-(34) 
for tan/? = 30 and 7 = 0, with the same convention in Table 1. 

behaviour of R 2b for large tan (3 follows from the complementary behaviour 
of its 0i and sin /3(pi + cos/3y?2 compositions shown in Fig. 3 (right panel). 
It is with the a b = sin (3 factor in eq. (35) that such a compensation between 
its opposite CP-components occur. The variation of R 2c follows from the 
02 composition of H 2 Fig. 3 following the same lines of reasoning used in 
discussing R lb above. 

Finally, Fig. 8 shows the variation of R 3b (left panel) and R 3c (right 
panel) with 7 for tan/3 = 4 (solid curve) and tan/3 = 30 (dashed curve). 
For tan/3 = 4, both R 3b and R 3c remain around unity because of the fact 
that there is little H — A mixing and (7^ 33 ) 2 (0) = 1. That R 3b remains flat 
for tan/3 = 30 follows from the interplay between its CP=+1 and CP=-1 
components as in R 2b . On the other hand, R 3c , for tan/3 = 4 follows simply 
from its vanishing 2 and 100% sin/3</?i + cos/3y? 2 components in Fig. 4 (left 
panel). Since 2 composition is negligibly small for large tan (3 (right panel of 
Fig. 4), variation of R 3c is mainly dictated by its CP=-1 component (short- 
dotted curve in the right panel of Fig. 4). For example, R 2c = 1 at 7 = n 
just due to its 100% composition in Fig. 4. 

From the study of the decay rates of the Higgs particles to bb and cc pairs 
one concludes that their CP-impurity causes significant changes compared 
to the CP-invariant limit. In particular, one notes the enhancement in R 3c 
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Figure 7: i?26 (left panel) and i?2c (right panel) defined in (34) as a function of 7 
for tan (3 = 4 (solid curve) and tan (3 = 30 (dashed curve) . 

near 7 = point, which is one order of magnitude above its value in the 
CP-conserving limit. 



4 Conclusion and Discussions 

In this work we have studied the possible effects of the supersymmetric CP- 
violating phases on the neutral Higgs scalars of the MSSM. We have adopted 
effective potential approximation in computing the radiative corrections and 
have taken into account only the dominant top quark and top squark loops. 
We now itemize the main results of the work: 

• Radiative corrections induce an unremovable relative phase between the 
two Higgs doublets. This relative phase remains non-vanishing as long 
as the supersymmetric CP-violating phases are finite, and determines 



the dynamics of the electroweak phase transition JT9] 



The CP=+1 and CP=-1 components of the Higgs doublets are mixed 
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Figure 8: R^b (left panel) and R^ c (right panel) denned in (34) as a function of 7 
for tan (3 = 4 (solid curve) and tan (3 = 30 (dashed curve) . 

up due to the mixing terms in the scalar mass-squared matrix which 

1. are proportional to the relative phase between the two doublets, 

2. increase with increasing tan (5. 

• The lightest Higgs scalar remains essentially CP-even irrespective of 
the super symmetric CP phases. Therefore, the MSSM has a light 
CP=+1 scalar as in the CP-respecting case which can, however, be 
distinguished from that of the CP-invariant theory, for example, by its 
reduced decay rate to bb pairs. 

• As in the CP-invariant case, there are two heavy scalars which have 

1. definite CP quantum numbers for small tan/9 values, 

2. no definite CP characteristics for large tan/5. 

• The strong mixing between the heavy scalars affect significantly their 
decay rates to fermion pairs, and such mixings can be important in 
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other collision processes testable at future colliders. Especially the 
decay rate of the would-be CP-odd scalar gets enhanced for < 7 < 
tt/2. 

• The supersymmetric CP-violating phases not only cause the creation 
of CP-violating mixings but also modify the couplings and masses com- 
pared to ones in the CP-invariant case. 

In the light of these results one concludes that the supersymmetric phases 
can be useful tools for obtaining manifestations of supersymmetry through 
their effects on collision and decay processes testable in near future colliders. 
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